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Based on the perturbative renormalization group (PRG) approach, we have studied 
dimensional crossovers in Hubbard ladders coupled via weak interladder one-particle 
hopping, t±. We found that the one-particle crossover is strongly suppressed through 
growth of the intraladder scattering processes which lead the isolated Hubbard ladder 
system toward the spin gap metal (SGM) phase. Consequently when t± sets in, there 
exists, for any finite intraladder Hubbard repulsion, U > 0, the region where the two- 
particle crossover dominates the one-particle crossover and consequently the d-wave 
superconducting transition, which is regarded as a bipolaron condensation, occurs. By 
solving the scaling equations for the interladder one-particle and two-particle hopping 
amplitudes, we give phase diagrams of the system with respect to U, t± (initial value 
of t±) and the temperature, T. We compared the above dimensional crossovers with 
those in a weakly coupled chain system, clarifying the difference between them. 

KEYWORDS: doped ladder, dimensional crossover, perturbative renormalization-group, spin gap metal, 
bipolaron condensation, d-wave superconductivity 



§1. Introduction 

Magnetic and electronic properties of ladder materials have attracted great interest. |l 
Central to these issues are the effects of the unusual spin-liquid state with a spin excitation 
gap in the undoped parent system on the electronic conduction in the doped system. Last 
year Uehara et al. j2|] discovered a superconductivity signal in the doped ladder system, 
Sr 14 _ x Ca a; Cu2404i, under pressure. The compound consists of the alternating stacks of planes 
with the CUO2 chains and with the CU2O3 two-leg ladders. The nominal valence of Cu is 

* E-mail:kishine@ims. ac.jp 



1 



+2.25, independent of x, so the system is inherently doped with 6 holes in the unit, which are 
mainly on the chains at x = 0. Increase of the low-energy spectral weight |J with increasing 
x indicates that the hole carriers are progressively redistributed from chains to ladders with 
Ca substitution for Sr. This feature was also confirmed theoretically through the ionic and 
cluster model calculations. [|J The behavior of the dc resistivity of the compounds changes 
from semiconducting-like to metal-like with Ca substitution, and superconductivity sets 
in for Sr .4Cai3.6Cu24O4i below T c = 12 K under a pressure of 3 GPa. || 

The spin excitation gap at the ladder cite in Sr 14 Cu240 4 i has been observed through 
NMR shifts and rates, |||7| and neutron scattering experiments. [|J A remarkable feature is 
that the spin gap in the ladder survives upon Ca substitution, |],[TIJ indicating the doped 
ladder compounds, Sr 14 _, r Ca. r Cu2404i, are in the spin gap metal phase which reflects the 
one-dimensional character of the ladder. According to the electronic structure calculation of 
Sr 14 _ a .Ca a; Cu2404i under ambient pressure within the local-density approximation, ]nj] the 
interladder hoppings are 5-20 % of the intraladder ones, indicating a pseudo-one-dimensional 
character of the system. 

The superconducting transition under pressure suggests that interladder one-particle hop- 
ping enhanced by applied pressure plays an important role. Recent experiments of the 
resistivity along the ladders, p c , of the single crystal S^.sCan.sCu^C^i by Akimitsu et 
al. [ 12 1 shows that the superconductivity sets in below 10 K under 3.5 GPa 
the temperature dependence of p c changing gradually from T-linear to T 2 



8 GPa with 
The ratio of the 

normal state resistivity, p c j p a (p a is resistivity perpendicular to the ladder direction in the 
ladder plane) also indicates a dimensional crossover from ID to 2D with increasing applied 
pressure. (12| 

The low-energy asymptotics of the isolated ladder system has been extensively studied 
mainly based on the renormalizat ion-group approach, with and without the aid of the 
bosonization technique, to the two-leg Hubbard ladder [|T^|T^,|T3|.|TB[[T^.|TE|,|T^] and the two-leg 
t-J ladder. E0],|2T|,|22|,E3l In the isolated Hubbard ladder, the most relevant phase is char- 



acterized by a strong coupling fixed point and is denoted by "phase I" by Fabrizio |TB[ and 
"C1S0 phase" by Balents and Fisher. [|16[] In this phase, only the total charge mode remains 
gapless and consequently the (i-wave superconducting correlation becomes the most domi- 
nant one [p^,[l5],[T6],|l7|] and the 4fcp-CDW correlation becomes sub-dominant one fl5l, |i6i|T7l 
at least when the intraladder correlation is weak. From now on, we call this strong coupling 
phase a "spin gap metal (SGM) phase". Within the bosonization scheme, exponents of the 
ci-wave superconducting and the 4/cpCDW correlations, K$ and Kc, satisfy the duality rela- 
tion, Ks ■ Kc = 1, suggesting the low-energy dynamics in the SGM phase is described with 
"bipolarons" , each of which lies along the rung of the ladder. |21],|24| The enhancement of 
the d-wave superconducting channel has also been found by numerical studies on the two-leg 
Hubbard ladder through the density matrix renormalization group (DMRG) analysis, |25|] 
the exact diagonalization study, |26| and the Monte-Carlo simulation. pTI] The exact diag- 
onalization study on the two-leg t-J ladder also gave evidences of a gapless total charge 
mode [28] and a <i-wave RVB state [2i,^9| in the doped system. 

So far effects of the weak interladder hopping have been studied through mean field approx- 
imations |J0],|3T| and power counting arguments. |32j Recently the present authors discussed 
effects of the interladder one-particle hopping, t±, on the low-energy asymptotics of the 
weakly coupled Hubbard ladder system, based on the perturbative renormalization-group 
(PRG) approach and discussed dimensional crossovers in the system. |B3 
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Although the dimensional crossover problem in the coupled ladder system is rather new 
aspect, the similar problem in the coupled chain system has a history of about twenty years. 
[33],|35|,|36| Effects of strong intrachain quantum fluctuations on the competition between one- 



particle and two-particle crossovers were first addressed by Brazovskii and Yakovenko. p7 



It was pointed out by Suzumura, |38|] based on the bosonization technique, that the two- 
particle crossover dominates the one-particle crossover only in the presence of large intrachain 
correlation. Recently the competition was again discussed in terms of the stiffness of the 
Tomonaga-Luttinger liquid. [BU 



Bourbonnais and Caron, on the other hand, first discussed the problem based on the PRG 
approach, [iD,^TL|2| where the intrachain interaction and interchain one-particle hopping 



are treated as perturbations to one-dimensional free fermion gas. In their formulation, the 
two-particle processes are generated by the interchain one-particle hopping in the course of 
renormalization. Based on their formulation, the dimensional crossovers in the weakly cou- 
pled chains were discussed in terms of the anomalous dimension of the Tomonaga-Luttinger 
liquid. || 

In the present paper we extend discussions on the results we shortly presented in our 
previous paper. |33j Based on the PRG approach, we treated the intraladder interaction and 



the interladder one-particle hopping as perturbations to the free fermion gas on the ladder 
and obtained the phase diagram of the system. In this extended version, a full account of 
application of the PRG scheme to the dimensional crossover problem in the weakly coupled 
ladder system is given. Furthermore, we re-examine the dimensional crossovers in the weakly 
coupled chain system, clarifying the difference between the coupled ladders and chains. 
We show that the main difference comes from the different behaviors of the scaling flows 
of respective coupling strengths, i.e, the different universality classes of the corresponding 
isolated systems. 

The outline of the present paper is as follows. In §2, we give a qualitative description 
of the notion of the one-particle and two-particle crossovers in the weakly coupled ladder 
system. In §3, we give a full account of formulation; the action, the scaling equations for the 
intraladder processes, the interladder one-particle process and the interladder two-particle 
process. In §4, based on the solutions of the scaling equations, we give phase diagrams of 
the system. In §5, we re-examine the dimensional crossover problem in the weakly coupled 
chains. In the final section, we give concluding remarks. Some technical details are left to 
Appendices. 

§2. Dimensional Crossovers: One-Particle vs. Two-Particle Crossovers 

The central problem here is how the weak interladder one-particle hopping, t±, affects the 
low-energy asymptotics of the system. We switch on the intraladder interaction and the 
interladder one-particle hopping, t±, as perturbations to the system specified by the intral- 
adder longitudinal (transverse) hopping, t(t') (see Fig. 1). |0[ Then two kinds of interladder 
processes, a one-particle process and a two-particle process, occur. We illustrate these pro- 
cesses in Fig. 2. The one-particle and two-particle processes correspond to the interladder 
hopping of the single-particle excitation and that of the two-particle (particle-particle or 
particle-hole pair) excitation, respectively. As the temperature scale decreases, dimensional 



crossovers are induced by a one-particle process or a two-particle process. |3^,|38],|39],[42| If 
the former dominates the latter, an interladder coherent band motion occurs. Then the 
character of the ladder in one-dimension is completely lost. On the other hand, if the lat- 
ter dominates the former, interladder coherent propagation of a pair of composite particles 
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occurs and the system transits to a long-range-ordered phase in the corresponding channel. 
In Fig. 2, we show the two-particle hopping in the d-wave superconducting channel, which 
corresponds to the interladder Josephson tunneling of bipolarons. We will see this is the 
most dominant two-particle process. 

The isolated doped ladder system has a spin excitation gap, which corresponds to the 
binding energy of the bipolaron. Thus the amplitude of the interladder tunneling of a bipo- 
laron is roughly estimated as t 2 ± /A a , where A CT is the characteristic energy scale of the spin 
gap. On the other hand, the amplitude of the one-particle hopping is t±. Consequently, the 
competition is determined by the ratio, t±/A (T , which is strongly affected by one- dimensional 
quantum fluctuations in the system. 

§3. Scaling Equations for Interladder One-Particle and Two-Particle Hopping 
Amplitudes 

To study the competition between the one-particle and two-particle crossovers, we set up 
scaling equations for the interladder one-particle and two-particle hopping amplitudes. 

3.1 Action 

We start with the path integral representation of the partition function of the system, 
Z = j Ve s , where the action consists of four parts, 

S = + + + S[ 2) , (3.1) 

with Sj 2) , andSjH eing the actions for the intraladder one-particle hopping, intral- 
adder two-particle scattering, interladder one-particle and interladder two-particle hoppings, 
respectively. T> symbolizes the measure of the path integral over the fermionic Grassmann 
variables. 

The intraladder one-particle process is diagonalized with respect to the bonding (B) and 
antibonding (A) bands. As shown in Fig. 3(a), we linearize the dispersion along the legs on 
the bonding and antibonding Fermi points, ±fcp m (m = A, B). In Fig. 3, R and L denote 
the right-moving and left-moving branches, respectively. The actions for the intra- and 
interladder one-particle hopping processes are written as 

5 ii 1) = E E E [^(^ii)CW^(iO + 5^n)c W^W] . (s- 2 ) 

K m=A,B cr 

and 

4 1} = "E E T, £ ±rn(k ± )[Ll a (K)LUK) + R* ma (K)RUK)}, (3.3) 
K m=A,B cr 

where R ma (K) and L ma (K) are the Grassmann variables representing the right- and left- 
moving electrons on the band m. The dispersions in the interladder action are given by 
£±A(k±) = +t±cosk± and £_lb(£;_l) = — t±cosk±. We use the notation K = (k\\,k±,ie n ) 
and K\\ = (fc||,i£ n ) where e n = (2n + l)ir/(3 is a fermion thermal frequency. We denote the 
momenta along the leg and the rung by k\\ and k±, respectively. The intraladder Green 
functions are written as 

g um (K\\) = [ie n -e um (hT\ (3.4) 
where the linearized dispersions are given by eR m (k\\) = VF{k\\—kF m ) and £L m (A;||) = vp{—k\\ — 
kprn)- The Fermi momenta of the bonding and antibonding bands are given by Ufb = kp + 
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t'/vp and k F A = kp—t'/vp, respectively. The Fermi velocities in principle depend on the band 
index as Vp m = 2t sin kp m , but we assume throughout this work that Vp m = vp and drop the 
band index, since the difference in the Fermi velocities does not affect the asymptotic nature 
of the SGM phase at least for small t'/t. |13],[16| In all the four branches (LB, LA, RA, RB) 
of linearized bands, the energy variables, e vm (v = R,L;m = A,B), run over the region, 
—E/2 < e vm < E/2, with E denoting the bandwidth cutoff. 

The intraladder Hubbard repulsion generates scattering processes depicted in Fig. 3(b). 
The action for the intraladder two-particle scattering processes is written as 

(2) _ 2lXVp \ - 01020-304 f* p* p j I \ - 01020-304 r* p* p_ j 

°\\ a 2-^i yo mo"i mo"2 m(7 3 ^04 ~T~ q / , yf ^mai^fh^ mcr 3 mcr ± 

I ^TTVp \ " O! 020304 T* p* p_ T f) 

where the summations are taken over band indices, spins, intraladder momenta and frequen- 
cies as 

E ^mi<7i -^m 2 ff2 -^m3ff3-^m 4 ff4 

= E E E S(Ki + K2-K 3 - ^4)^^(^)^(^2)^3.3(^3)^404(^4) 

Ki,— ,K A mi,-,m4 01, ■■■,04 

with m and fa being different band indices. The dimensionless intraladder scattering strengths 
are given by 

0-JO20-3O4 _ r _ (2)r Sj /o a\ 

<J[i y^t u aia3 u a2(74 i/^ u oi 0-4^2 0-3 ; V°- u / 

where dimensionless quantities gjp and denote backward and forward scattering strengths, 



respectively, with the flavor indices, [13] fi = 0, f,t. The usual scattering strengths with di- 



mension of the interaction energy are 2i\vpg^ L ' . We neglect the interband backward scattering 



terms like 9b 1<J2 ' Tia4 L*Bo 1 R*A C j 2 RBa 3 LAa i on the ground that these processes do not seriously 



modify the asymptotic nature of the C1S0 phase for finite t' (13|,[T6]]. Since we neglect both 
of the band index dependence of the Fermi velocities and the interband backward processes, 
the intraladder transverse hopping, t' , never appears explicitly in the present work. 

The action for the interladder two-particle hopping processes are decomposed into CDW/i, 
SDW/U, SSfi (singlet superconducting) and TSfi (triplet superconducting) channels (fi = 
0, /, t specify the corresponding flavor indices) as 



o(2) ~F \ \ * [T/DW/omm*//imm i T/DW/omirw/omm _j_ i/DW/nrann/'/i 

A DW=CDW,SDW 



?TVf x " }T [V^0£ m *0^ m + V t s O£ m *Og™ + VfOf m *Of m ~\ , (3.7) 



S=SS,TS 



where V^f denotes the interladder pair tunneling amplitude for channel M/z. m and m denote 
different band indices. The summations are taken over band indices, intra- and inter-ladder 
momenta and boson thermal frequencies, Q = (q\\, q±, iooi), with uii = 2lir//3, 

j2v^oz im2 *oz 3m4 = Y, E v^ M oz im2 *(Q)oz 3m4 (Q). 

Q mi,-,m 4 
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The corresponding composite field variables are defined by 

0^ w (Q)=p- 1/2 T,R*mAK + Q)L m , a (K), 

d^ w (Q)=(3- 1 / 2 Y,R™AK + Q)a^L m , a ,(K), (3.8) 
0^\Q) = f3~ 1 ' 2 J2 oR m(J {-K + Q)L m , 5 (K), 
6%f{Q) = /r 1/2 E °R™(-k + Q)^a'L m ^(K), 

where a are Pauli matrices, a = —o~, and the summations are taken over spins and K. 

3.2 Scaling of the intraladder two-particle scattering processes 

The idea of scaling is to eliminate the short- wavelength degrees of freedom to relate effective 
actions at successive energy scales. Based on the bandwidth cutoff regularization scheme, 
we parametrize the cutoff as E{1) = Eqc~ 1 with the scaling parameter, /, and study how the 
action is renormalized as I goes from zero to infinity. Details on the derivation of the scaling 
equations are left to the Appendices. 

Scaling equations for the intraladder scattering vertices, g^\ are generally written as [34] 

^|^-£(0-^,W. (3-9) 

Within the 3rd order PRG scheme, the first term on the r.h.s, Zy(l), comes from the vertex 
correction diagrams represented in Figs. 4(a) and 4(b) for i = 1 and 2, respectively. The 
field rescaling procedure (see Appendix A) requires the second term, 2z\\, to be 

*||(0 = 0b W (O a + flS°(0 2 + 9?\lT + 9?{l? + 9 { '\lf + 9?\l? 

-9 { o\i)g?\i) - gf\i)g?\i) - gPflgPO), (3.io) 

which comes from the intraladder self-energy diagrams represented in Fig. 4(c). We give 
a full detail on the derivation of z\\ in Appendix B. The full expressions for the scaling 
equations, (3.9), are given by 

|«tf> (0 = -%? ,2 + M + 2( 9 <v > 2 - s < W) 



(1)/ (l) 2 . (I) 2 . (I) 2 . (2) 2 (1) (2)s /oii\ 

-2^o (9o +g) +9t + 9t ~ 9t 9t h I 3 - 11 ) 



d 



2)/ n / l) z . (1 Y . (2 z x . Q / 2) (2) z 1) (2 (2 K 

(0 = -m +gl +9t ) + 2 (g) 9t - m 9t g) ) 

+2&*gf - g^ - gfg^ - g? 2 gf - 2gf\g^ + gf - g? g?\ (3.12) 
9?\l) = -2( 9 f + 9 U 2 -gM 9 V) 

. r>( (1) (2) 2 (1) (1) (2)s Q (1)/ (l) 2 . (I) 2 . (I) 2 . (2) 2 (1) (2)x /o i o\ 

+2(ffo ft ~% 9t 9t )-^9) \9q +9f + 9t + 9t ~ 9t 9t ), (3-13) 

dV 



dl 



±g?(D = -{gf - gf )2 ) + 2(9^9^ - gl 1] gl 2) gi 2) ) 

2 to\ m2 ci\ n\2 m n\3 _ in\ , n\2 fo\2 



l0 (l) z (2) (1)^ (1) (1)^ (1) (l)- 3 (2)/ (l) z . (2) z (1) (2)x foiA\ 

+m % - 9% % - % 9) - 9) - 2 g) m +gl - m gl '), (3.14) 

dl 



= -2(2g?gf - g?gf - gfg? + + ^ 
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+2{2^gVgf-$<Pgf-gytf>ti 



(l)«(l)«(2)^ 



(3.15) 



,(2)' 



d (2) 



eft 



9f 



{2y 



9i 



9i 



(2)< 



(1) (2) (1) (2) (1) (2) > 

#o £o ~9f 9f ~ 9t 9t ) 



^U 2 ) (2) 



^U 2 ) (2) 



-2gf\g^ + £>- + + gf~ + ^~ + g\ 



,(2) 



(1)' 



,(2)' 



(1) (2) (1) (2) (1) (2)n 

£o #o - 9/ 9/ ~ 9t 91 , 



(3.16) 



where we omitted the argument I in the scattering strengths on the r.h.s. The same equations 
are obtained by setting = g^ = in Eq.(A5) of Ref.[13]. 
Starting with the Hubbard type initial condition, 



g^(0) = U=U/47rv F >0, 



the scaling equations lead to the strong coupling fixed point (13 

0. 



(i)* 

(2)* 

9o = 



-1, 

3-2U 
4 ' 



(1)* 

9f 

_(2)* _ 1+2U 

9f ~ 



9i 



(i)* 



9i 



(2)* 



1- 
1. 



(3.17) 



(3.18) 



To obtain a physical description of the fixed point, it is useful to analyze the isolated 
Hubbard ladder system by means of the bosonization technique. fTS], [Tj], ITS], |T6l, [T7|, [T8|, |3"8"| . |2T)| . 

^-v sf\ I ,-v» \ r\ vi s~\ *• s~i t-1 I rm \ c\ -v r-\ 



21], ^2|, Q We introduce the phase fields, 0^(x) and 6 u (x), where <j) v {x) and -d x 9 u (x) are 



conjugate, 



7T 



(3.19) 



The mode-indices, /x, z/ = pA, pi?, crA, cri?, represent the charge- and spin-modes on the 
bands A or 5, respectively. The charge- and spin- modes are combined into the total-mode 
(0-mode) and the relative-mode (7r-mode) as 



■>pB 



± 



tB 



± 



(3.20) 



where + and — denote the 0-mode and 7r-mode, respectively. Then the "non-interacting " 
part [including the g^ processes (/i = 0, /) for i = 1 with parallel spins and those for i = 2] 



of the Hamiltionian for the isolated Hubbard ladder system is written as fT3|,|T5|. p"6 

Ho - 



£ It 



2tt 



u v K v (d x 6 u ) 2 + (d^) 2 
Ay 



The velocity, u u , and the stiffness, K u , of each mode are given by 



(3.21) 



V F \ 1 



9l, K » 



< l + 9u 
1 - 9v 



where 



g p± = (gW - 2^) ± (g) 



o„(2)\ - _ i „(1) 

^5/ J, fl'o-i - 9o ± 9f ■ 



(3.22) 



(3.23) 
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The low-energy asymptotics of the isolated Hubbard ladder system is characterized by the 
jrong couplin 
correspond to 



strong coupling fixed point of three scattering strengths, gl = g^* = —g^* = 1, which 



= K a+ = K a — = 0. (3.24) 



The scale-invariance of K p+ — y(l — 2U)/(1 + 2U) means that the total-charge mode re- 
mains gapless. On the other hand, 1/K p - = and K a+ = K a ^ = mean that the phases 8 P - 
and 4> a ± are respectively locked, suggesting the corresponding modes acquire a gap. [|l~3l.|T5| 
Therefore, in this phase, only the total charge mode remains gapless and consequently the 
(i-wave superconducting correlation becomes the most dominant one [|1^, [15], [16], [17]] and the 
4/cp-CDW correlation becomes sub-dominant one [15, |16|,|T7|1 at least when the intraladder 



correlation is weak. The phase characterized by this fixed point is denoted by the "phase I" 



by Fabrizio |13[ and the "C1S0 phase" by Balents and Fisher. JTfj] Throughout the present 
paper we call this phase a spin gap metal (SGM) phase. Exponents of the o?-wave supercon- 
ducting and the 4fcp-CDW correlations, Ks and Kc, in the SGM phase satisfy the duality 
relation, Ks ■ Kc = 1, suggesting the low-energy dynamics in the SGM phase is described 
with bipolarons, each of which lies along the rung of the ladder. pT].[24[ 



3. 3 Scaling of the interladder one-particle hopping amplitude 

The scaling of the interladder one-particle process is fully determined by the intraladder 
self-energy effects. The scaling equation, represented in Fig. 5, is written as 

—£i =!-,„(/). (3.25) 
We solve the equation with the initial condition, 

t±(0)/E = t ± o- (3.26) 
The scaling equation, (3.25), and the fixed point, (3.18), lead to 

d1nt±(l) z^oo 



dl 



-U 2 /2 - 7/8, (3.27) 



and consequently t±(l) becomes always irrelevant at the final stage of the scaling procedure. 
However, at an early stage of scaling, the r.h.s of eq.(3.25) remains positive since the intral- 
adder couplings do not grow sufficiently as yet and consequently t±(l) grows. The r.h.s of 
(3.25) changes its sign from positive to negative in the course of renormalization and then 
t±(l) begins to decrease and is finally scaled to zero. 

If, in the course of the scaling, tj_(l) = t±(l)/E attains an order of unity around some 
crossover value of the scaling parameter, Z cross , qualitatively specified by 

t± (/cross) = 1) (3.28) 

the weakly coupled ladder picture breaks down. Then, the notion of "relevance" or "ir- 
relevance" of the interladder one-particle hopping loses its literal meaning, since the PRG 
treatment for t± becomes inapplicable for I > Z cr0 ss- If Across precedes the scaling parameter 
specifying the two-particle crossover, the interladder one-particle hopping becomes coherent 
and the system is regarded as scaled to a "two-dimensional" system via the one-particle 
crossover. p0[|41|,p 



In Fig. 6(a) we illustrate the scaling flows of t±(l) with a fixed value of the intraladder 
repulsion, U = 0.3, and t± = 0.01, 0.02, 0.03, 0.04. The peak positions are independent 
of t± and the peak heights are proportional to t± . In Fig. 6(b) we illustrate the scaling 
flows of t±(l) with a fixed initial value of the interladder one-particle hopping t±o — 0.01 
and U = 0.1,0.2,0.3,0.4,0.5. The peak positions and heights are sensitive to U. We see 
from Fig.6 that for small U, t±(l) reaches unity even for small t±o, while for larger U, t±(l) 
are strongly suppressed and very large t± is required for t±(l) to attain an order of unity. 
As the intraladder correlations become stronger, the one-particle process is more severely 
suppressed. 



3.4 Scaling of the interladder two-particle hopping amplitudes 

The scaling equation for the interladder two-particle hopping amplitudes, , are gener- 
ally written as 



dl 



/"(0 + *"(0O). 



(3.29) 



where fjf denote the generators of V™ for channels Mfi (M = CDW, SDW, SS, TS and 
[A — 0, /, t). Note that zf^ has [V^ M ] n contributions with n = and n — 1. In Figs. 7(a) 
and 7(b), we show the diagrammatic representation of fjf for M = CDW/SDW and SS/TS, 
respectively. In Appendix C, we give the details on the derivation of fff . 

The leading-order scaling equations for V^f are diagrammatically given in Figs. 8(a) and 
8(b), for M =CDW/SDW and SS/TS, respectively. Full expressions for the scaling equa- 
tions, ( ET29D , are written as 



dV DW (l) 
dl 

dvp w {i) 

dl 



dl 

dl 
dVf(l) 



\ [t ± (/), DW (/)] 2 cosg ± + So DW (0^o DW (0 - 



1 ~ 



-?±(o 2 K w (0 2 +^ w (o 



DW/ 



cosgj_ 



+ g^(l)Vf(l) + g^ w (l)V f uw (l) - - [V t DW (lj 2 + V f uw (l) 



rDW/ 



1~ 



i-t ± (l)g? w (l)gf w (I) cos g± 



DW/ 



(0 + #"(0^(0 - vr w (0^/(0, 

-\h{lf [gM + <7 t S (0 2 ] cosg ± + g s (l)V s (l) + rf(0*?(0 

' \ s (i) 2 + v t s (if 



dl 

4 s 



-L {if g%lf cos qi _ + g S f (l)V f s (l) - \vf{l)\ 



(3.30) 
(3.31) 

(3.32) 
(3.33) 

(3.34) 



= -lt ± (lT9 S o(l)9?(l) cosq ± + g*{l)V*{t) + gfV s (l) - V S (l)V t s (I) , (3.35) 



with DW=CDW/SDW and S=SS/TS. Coupling strengths of the composite particles in these 
channels are given by 



^(0 = ^(0-2^(0, 



.SDW/ 



^•'(0 = »i a) (0, 



gf(i) 



(3.36) 
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We here introduce the composite fields of the s- and d-wave superconducting channels by 
T|,|lJ,|l6l[Ij],|23,fg,H 



SCs 



mBB i mAA /n mBB mAA 



'SS 



ss • 



(3.37) 



Cooper pairs in the d-wave channel (=bipolaron on the rung) are formed on the adjacent legs on 
the same rung (see Fig. 2). The corresponding two-particle hopping amplitudes are constructed 

as 



V SCs (l) = - V^(l) + Vf a (Z) , V^ d (l) = - V^(l) - V™(1) 



1 r 



rSS/ 



rSS/ 



SCd/ 



1 r 



rSS/ 



that satisfy the scaling equations 

d ySCs 



dl 
d ySCd 



dl 

with the coupling strengths 



t ± (l)9 SCs (l)fcos q± + 2gl c °(l)V SSs (l) - \v SSs (ir, 
t ± (l)9 SCd (l)} 2 cos q ± + 2gl Cd (l)V SCd (l) - \v SCd (l)\ 



1 r 



which are scaled to the fixed point, 

,sc s * = _(i + 2U) /8, ^ SCd * = (15 - 2U) /8. 
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We have solved simultaneously the scaling equations, (3.11) 
(3.39) and (3.40) with the initial conditions 



(3.16), (3.25), (3.30) 



g^(0) = U/47rv F = U, t ± {Q)/E 



t±o, 



C(o) 



0. 



(3.38) 

(3.39) 
(3.40) 

(3.41) 

(3.42) 
(3.35), 

(3.43) 



We put q± = and q± = n for the superconducting and density-wave channels, respectively. In 
Figs. 9(a) ~ (c), we show the scaling flows of the coupling strengths for the composite fields, (3.36), 
and V^ 1 (/) for U = 0.2,0.3,0.4 and t±Q = 0.01, where the vertical broken line corresponds to 
the scaling parameter at which V SCd diverges. We clearly see that the divergence of V always 
occurs at the highest energy scale. This situation is quite reasonable on physical grounds that 
the interladder pair tunneling stabilizes the most dominant intraladder correlation, i.e, the d-wave 
superconducting correlation. As is seen from Figs. 9(d) ~ (f), the competition among the channels 
other than the SCd channel is rather subtle. Below we focus on the d-wave superconducting channel. 

The first, second and the third terms on the r.h.s of eq.(3.40) play separate roles in the early, 
intermediate and final stages of the scaling of V SCd . At the early stage, the first term generates 
a finite magnitude of V SCd . At the intermediate stage, the second term induces an exponential 
growth of V SCd . At the final stage, the third term causes divergence of \/ SCd at a critical scaling 
parameter l c defined by 



V SCd (l c ) 



-oo. 



(3.44) 



§4. Phase Diagram of Weakly Coupled Hubbard Ladders 

Based on the scaling flows of t±(l) and V SCd , we show the competition between the one- and 
two-particle crossovers. The scaling flows of the intraladder system toward the SGM phase are best 
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visualized through the scaling flow of K a+ , the stiffness of the total-spin mode. We here again note 
that K a+ = corresponds to the fully developed spin gap at the low-energy limit. In Fig. 10, we 
illustrate the scaling flows of K a+ , t±(l) and V SCd with various initial conditions. We see from 
Fig. 10 that, for (U,t± ) =(0.3,0.01), (0.3,0.02), (0.4,0.01) and (0.5,0.01), the two-particle crossover 
occurs at I = l c . On the other hand, for (U,t±o) =(0.3,0.03), (0.3,0.04), (0.1,0.01) and (0.2,0.01), 
t±(l) exceeds unity at l cmss before V SCd diverges. Then the one-particle crossover occurs. 

The scaling parameter is identified with the absolute temperature as Z = ln-^n. Thus we de- 
fine the one-particle crossover temperature, T cross , and the ci- wave superconducting transition 
temperature, T c , as 

r T 1 — Z7 1 ci Across 

cross -MJ C i {A -\\ 

T c = E e- l «. 

We obtain phase diagrams of the system with respect to (U, tj_o) and the reduced temperature, 
f = T/Eq. 

4-1 t±o-T phase diagram 

First we show, in Fig. 11, the phase diagram spanned by t±o and T for U = 0.3. Roughly 
speaking, we may regard increasing tj_o as applying the pressure, under which the bulk supercon- 



ductivity was actually observed in a doped spin ladder, Sri4_ x Ca IE Cu2404i. [g, 12| We found that 
there exists a crossover value of the interladder one-particle hopping, tj_ c ~ 0.025. 

For < tj_o < t± c , the phase transition into the ci-wave superconducting (SCd) phase occurs at 
a finite transition temperature, T c , via the condensation of bipolarons driven by the two-particle 
crossover. In the temperature region, T < T c , interladder coherent Josephson tunneling of the 
bipolarons occurs. Now we must use caution in identifying the finite temperature phase above T c , 
where the system is in the isolated ladder regime. As the temperature scale decreases, the isolated 
ladder systems are gradually scaled to their low-energy asymptotics, the SGM phase, toward the 
zero temperature. The gradual change of darkness in the SGM phase in Fig. 11 schematically 
illustrates this situation. The SGM phase is characterized by the strong coupling fixed point of 
intraladder scattering strengths, # t = gl = —g^ = 1, or l/K p _ = K a+ = = (see 

(3.18) and (3.24)). We see from Fig. 10 that the critical scaling parameter l c is always in the region 
around which K a+ almost reaches its fixed point value. Thus we expect that the spin gap is well 
developed near T c and the phase transition at T c can be identified with the transition from the 
SGM phase to the SCd phase. Within the framework of the PRG approach, however, we cannot 
say for certain whether the spin gap survives in the SCd phase or not. 

For t± c < t±Q, the system undergoes the crossover to the 2D phase via the one-particle crossover. 
The temperature scale of the one-particle crossover is not low enough to ensure a well-developed spin 
gap since the one-particle crossover takes place at an early stage of scaling, where the intraladder 
couplings are still far away from their strong coupling values, as can be seen from Fig. 10. Thus 
it is disputable to assign the phase above T cross to the SGM phase. In the temperature region, 
T < T croaa the transverse coherent band motion occurs. In section 4.4, we give brief comments on 
the physical nature of the 2D phase. 

In this paper, we consider the case where t is much larger than t±Q and is not so large as only one 
band to cut the Fermi surface (these conditions are actually satisfied in the real Sr 14 _ x Ca x Cu2404i 
compounds 0]). Then, the (hypothetical) isolated ladder system is always scaled to the SGM 
phase. When t' decreases and becomes comparable to t±o, the system may be regarded as a 
coupled chain system, since, in such a case, the alternation of the transverse hopping would be 



11 



unimportant. As we shall discuss in §5, in the coupled chain system, the one-particle crossover 
converts the system to the 2D phase as far as U is not extremely large. Then the SCd region would 
shrink with increasing t to be connected with the case of comparable t! and t±Q. However, our 
result (Fig. 11) would not be qualitatively changed. 

4-2 U-T phase diagram 

In Fig. 12, we show the phase diagram spanned by U and T for t±Q = 0.01. We see that there 
exists a crossover value of the intraladder repulsion, U c ~ 0.22. For < U < U c , the system 
undergoes the crossover to the 2D phase via the one-particle crossover, while for U c < U, the SGM 
phase transits to the SCd phase via the two-particle crossover. 

4-3 Crossover values of the interladder one-particle hopping: t± c 

In Fig. 13, we show how the crossover value, t± c , introduced in Fig. 11, depends on U. For 
£±o > tj_ c , the system undergoes the crossover to the 2D phase, while for t±o < t± c , the SGM 
phase transits to the SCd phase. t± c is always finite for U > and becomes zero only for U — 0. 
Thus we see that the crossover value, t± c , always exists for U > 0. The SCd phase vanishes 
only for {7 = 0, where the scaling equation, (3.25), gives T cross = t±Q. 

4-4 Remarks on the 2D phase 

Here we give a few remarks on the 2D phase. In the 2D phase the physical properties of the 
system would strongly depend on the shape of the 2D Fermi surface which may be characterized by 
the misfit of the rungs in the neighboring ladders which actually exists in real Sri4_ 2 ,Ca x Cu2404i 
compounds (see Fig. 14(a)). The band structure of the system is characterized by bonding (B) and 
antibonding (A) dispersions, 



/ k k 

e AB (k) = -2tcos£;|| ± Wt' 2 + 4t" 2 cos 2 y + 2ft" cos cos % (4.2) 

where the + and — signs are taken for A and B bands, respectively. It was pointed out by 
Yamaji [Effl that in the two-dimensional two-band system, a superconducting transition is possible 



via the so-called Suhl-Kondo mechanism, [48,49| where the interband exchange- like interaction is 
strongly enhanced by an interband nesting. Then the pairing interaction is caused by the processe 
expressed by the interband particle-hole ladder diagrams. Quite recently it was reported [^] that, 
by using a numerical fluctuation-exchange (FLEX) method, a superconducting transition actually 
becomes possible in a coupled Hubbard ladder system with the same configuration as that shown 
in Fig. 14(a), where the two-band structure of the system is taken into account. 

To check how the interband polarization is enhanced by the interband nesting, we calculate the 
static interband polarization function written as 

1 tanh[f(£ A (£+?/2) -//)]- tanh[f(e B (£-g/2)-A/)] ( . 
XA B (q\\,q±) = ~}2 2 -, , (4.3) 

2 j. e A {k + q/2) -£ B (k-q/2) 

In Fig. 14(b), we show Xab(c[\\, q±) on the 1st Brillouin zone along the line T(0, 0) — > X(7r, 0) — > 
M(7T, it) — ► r for for t = t' and t" jt = 0,0.2,0.3 with a chemical potential /i = and a 
temperature, (3~ l = T = 10~ 5 t. We see that an increasing t" decreases the degree of interband 
nesting and consequently the interband polarization function decreases. Then the exchange-like 
pairing interaction would also decrease. Thus, when a superconducting transition via Suhl-Kondo 
mechanism becomes possible in the 2D phase, the supercondicting transition temperature, T c , 
should decrease with an increasing interladder coupling. 
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§5. Comparison with Dimensional Crossovers in Weakly Coupled Chains 

In this section we compare the dimensional crossovers in the present system with those in the 
weakly coupled Hubbard chains. Although the PRG approach to the dimensional crossover problem 
for the weakly coupled chains has been extensively studied by Broubonnais and Caron, [42| it is 



instructive to re-examine their work, clarifying the difference in the nature of the dimensional 
crossovers in the two cases. 
The action for the weakly coupled chain system consists of four parts, 

C _ c(l) , C(2) , , C(2) (K i\ 

^chain — ^chain|| °chain|| ~T °chain_L °chain_L' \°- L ) 

where S^ainl] > ^chamH ' 'S'chainX an d 'S'chainX denote the actions for the intrachain one-particle hopping, 
intrachain two-particle scattering, interchain one-particle and interchain two-particle hoppings, 
respectively. As shown in Fig. 15(a), we linearize the dispersion along the chains on the two 
Fermi points, rtfejr. The intrachain Hubbard repulsion generates scattering processes depicted in 
Fig. 15(b). The processes are specified by dimensionless quantities g^ 1 ' and denoting backward 
and forward scattering strengths, respectively. [|34| The usual scattering strengths with dimension 
of the interaction energy are nvpg^. 

The action for the interchain two-particle hopping processes are decomposed into CDW, SDW, 
SS (singlet superconducting) and TS (triplet superconducting) channels and is written, instead of 
(3.7), as 

s$L± = -^ E E^d<?d-^ E E^s^s, (5.2) 

4 D=CDW,SDW Q 4 S=SS,TS Q 

where V M denotes the amplitude of the interchain two-particle hopping for channel M. The 
corresponding composite field variables are defined by 

Ccdw(Q) = P~ 1/2 E R t( K + Q)L a (K), 

O sdw (Q) = r 1/2 J2K(K + Q)cr aa ,L a ,(K), (5.3) 

Oss(Q) = p- 1/2 J2<rM-K + Q)La(K), 

Ots(Q) = /T 1/2 E <rR*(-K + Q)^L^(K), 

where a are Pauli matrices, a = —<J, and the summations are taken over spins and K with 
K = (fen, k±, ie n ) for each Q = (q\\, q± 7 iui) with fermion and boson thermal frequencies, £ n = 
(2n + l)ir/[3 and uji = 2ln/j3, respectively. We denote the momenta along and perpendicular to 
the chains by fey and k±, respectively. 

The diagrams which give the scaling equations for the intrachain scattering processes are the 
same as in Fig. 4, but some of them are canceled out to give, instead of (3.11) ~ (3.16), simple 



scaling equations [34| 



^ = -S (1 »(0 2 -^ (1) (f) 3 , (5.4) 

^ = -i 9 "W-i 9 (1 W. (5.5) 

Starting with the Hubbard type initial condition 

0«(O) = U/ttv f = AU > 0, (5.6) 
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the scaling equations lead to the fixed point 

g (D* = Q) g (2)* 



2U. 



(5.7) 



The fixed point characterizes the Tomonaga-Luttinger liquid (TL) phase. Since monotonically 
decreases, g^ 3 terms in (5.4) and (5.5) are not essential for the scaling property so that we can 
drop them. Then we obtain 



9 (1 \l) 



AU 



1 + AUl 



g (2) (l) = 2U + 



2U 



1 + AUl 



(5.8) 
(5.9) 



The TL phase is different from the SGM phase of the Hubbard ladder in that it belongs to 
the weak-coupling universality class and is gapless. Reflecting the weak-coupling nature, the 
intrachain self-energy effects are much weaker than in the ladder system. Consequently the strong 
suppression of the interchain one-particle process which is characteristic of the coupled ladders 
never occurs in the coupled chains. Furthermore, since the TL phase is gapless, the interchain 
two-particle hopping amplitude is ill defined (see the discussion given in §2). 

The scaling equation for the interchain one-particle hopping amplitude, t c hain_i_> [42] is, instead 
of (3.25), written as 



dlnt chlLin± (l) 
all 



^chain|| (J>) j 



with 



1 



g (l \i) 2 + g (2 Ki f - g (1 Kl)g (2 \l) 



•^chain|| (0 

Using (5.8) and (5.9), we obtain 

£chain±(0 = £chain±0exp 



■,(2) 



(1)/ 



,(2), 



[l-U 2 )l + 



3U 



4 Kl + AUl 



(5.10) 



(5.11) 



(5.12) 



Thus, contrary to the coupled-ladder case, t± becomes relevant for weak repulsion, U < 1, where 
the PRG scheme is reliable. We show this situation in the upper half plane of Fig. 16, where the 
scaling flows of £ c hain_i_(0 are shown for different U = 0.1, 0.2, 0.3 with t c hain±o — 0.01. The 
one-particle crossover temperature is defined by = e~' cross , where l^ss * s determined by 

^chain± Across 1 ) = 1- Hll ^ ne scaling behavior of t c h a in±(0 is also expressed in terms of the anomalous 
exponent, 6, as t ch ain±(0 = Win±o exp[(l - 9)1) which gives T c c r h ^ n = E [t cha , in±0 / £ ] 1/(1_e) - The 
exact solution 51, 3^,|53| tells us 6 satisfies 9 < 1/8, which again indicates t± is always relevant. 
Structures of the leading-order scaling equations for V M (M =CDW, SDW, SS, TS) are the 



same as those shown in Fig. 8, which give 42 



dV™(l) 
dl 



chain_L 



(% DW (0 



cosg ± + -g 



DW 



{l)V 



DW 



(I) - -V™{1)\ 



for DW=CDW or SDW and 
dV s (l) 



dl 



1 r~ 
t 



chain_L 



cosq ± + -g s (i)V s (i)--V s (lf 



(5.13) 



(5.14) 
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for S=SS or TS. We put q± = and q± = 7T for the superconducting and density-wave channels, 
respectively. Coupling strengths of the corresponding composite particles are given by 



g<*™(l)=gW(l)-2gP>(t), g SDW (7) = # (2) (/), 
» ffl (0 = -a/ (1) (0+ff (2, (0), ff TS (0 = ff (1) (0-» (a) (0- 

In the lower half plane of Fig. 16, we show the scaling flows of \/ SDW . 

In Figs. 17(a) and (b), we show the phase diagrams of the weakly coupled chains for U = 0.3 and 
for t c hain±o_ = 0.01, which correspond to Figs. 11 and 12 in the weakly coupled ladders, respectively. 
For small U and £cham_l_o> where the PRG scheme is reliable, the system always undergoes a crossover 
to a two-dimensional phase (2D phase) via the one-particle processes. We show in Fig. 17 for 
guidance the temperature, Tsdw> at which y SDW diverges. 

For the reasons already stated in §4.1, it needs care to identify the finite temperature phase above 
-^cross 1 - We can identify it with the TL phase characterized by the fixed point (5.7) only for very 
small £ c hain_i_o or large U, where the one-particle crossover occurs at a very low temperature scale. 
For large £ c hain_i_o or small U, the one-particle crossover occurs at a very early stage of the scaling 
where the intrachain scattering processes are still far away from their fixed point corresponding to 
the low-energy asymptotics of the system. The ambiguity in regarding the phase at T > T^f s n as 
the TL phase has some relevance to well known objections against the PRG approach to the present 
problem, as first emphasized by Anderson. |54|] He stressed that it is crucial to treat the intrachain 
interaction exactly before switching on the interchain one-particle hopping. The PRG scheme 
misses this point, since we there treat the interaction and the interchain one-particle hopping, on 
equal footing, as perturbations to the one-dimeniosnal free fermion system. Along Anderson's 
claim, by using the exact propagator of the TL liquid and calculating the momentum distribution 
function, Castellani et al. demonstrated that the TL phases are unstable with respect to 
arbitrary small interchain one-particle hopping, even if the interaction is strong. Their conclusion 
supports the results obtained through the PRG scheme, contraly to the expectation of Anderson. 

On the other hand, the effects of the interchain one-particle hopping have been studied in view of 



"coherence" or "incoherence" rather than its "relevance" or "irrelevance". [ 56 , 57 1 With the aid of 



the knowledge of the two- level system, it was demonstrated 56] that there exists a critical value of 
the inter-Luttinger-liquid hopping below which there is no coherent one-particle hopping between 
the liquids, which contradicts to the PRG results. At present, however, it remains an unsettled 
question whether the TL phase is stable against the interchain one-particle hopping or not. 

§6. Concluding Remarks 

In the present paper we applied the perturbative renormalization group (PRG) method to discuss 
dimensional crossovers in Hubbard ladders coupled via weak interladder one-particle hopping, t±. 
We set up and solved the scaling equations for the interladder one-particle and two-particle hopping 
amplitudes by treating the intraladder interactions and the interladder one-particle hopping as 
perturbations to the free electron system on the isolated ladders. We found that the scaling 
flow toward the spin gap metal phase (SGM phase) in the isolated Hubbard ladder strongly 
suppresses the interladder one-particle process and consequently, for any finite intraladder 
Hubbard repulsion, U > 0, there exits a finite crossover value of the interladder one-particle 
hopping, t± c . 

For < t± < t± c , the bulk d-wave superconducting phase is stabilized via the condensation 
of bipolarons, where interladder Josephson tunnleling of bipolarons occur. The superconducting 
transition occurs at the low temperature scales where the intraladder scattering processes nearly 
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attain their strong coupling values (see Fig. 10), which suggests the spin gap is already well developed 
around the temperature region. For t± c < t±, the system undergoes the crossover from the spin gap 
metal phase to the two-dimensional phase (2D phase) where the interladder one-particle hopping 
becomes coherent and one-dimensionality of the isolated ladder is completely lost. The temperature 
scales of the one-particle crossover is, however, not low enough to ensure a well-developed spin gap 
(see Fig. 10), since the one-particle crossover occurs at an early stage of scaling, where the intraladder 
system does not approach the strong coupling region as yet. Thus it is disputable whether we can 
assign the region T > T cross to the SGM phase or not. 

Our results are best summarized in the t±-T phase diagram, Fig.ll. Roughly speaking, we 
may regard increasing t± as applying pressure. In the doped spin ladder, Sr 2 .5Can.5Cu2404i, the 
superconductivity sets in below 10 K under 3.5 GPa ~ 8 GPa. The observed transition temperature, 
T° bs , increases with the pressure and reaches the maximum 6K around the optimal pressure P opt ~ 
4.5 GPa. [ 12 1 Our result in Fig.ll qualitatively reproduces the increase of T° bs for P < P Q pt- For 
P > -Popt, the T c obs gradually decreases and the resistivity along the ladder, p c , shows gradual 
change from T-linear to T 2 -dependence with increasing the applied pressure. [O] This fact strongly 



suggests that for P > P op t the system is regarded as an anisotropic 2D Fermi liquid. In our scheme, 
the existence of t± c would correspond to the existence of P pt- As was suggested in section 4.4, 
it is reasonable to suppose that in the 2D phase a superconducting instability via Suhl-Kondo 
mechanism j|7j] becomes possible due to the interband exchange-like interaction which is strongly 
enhanced by an interband nesting. Then the applied pressure tends to decrease the degree of the 
interband nesting and consequently the transition temperature should decrease with increasing the 
applied pressure. 

We also discussed difference between dimensional crossovers in weakly coupled Hubbard chains 
and ladders within the framework of the PRG. The difference originates from different universality 
classes to which the corresponding isolated systems belong. The isolated Hubbard chain belongs to 
the weak coupling Tomonaga-Luttinger liquid phase which is gapless, while the isolated Hubbard 
ladder belongs to the strong coupling phase with the spin gap. In the former case, the two-particle 
process is always dominated by the one-particle process, unless we assume very large intrachain 
Hubbard repulsion which is out of the perturbative scheme. On the contrary, in the latter case, 
the one-particle processes are strongly suppressed through growth of the intraladder scattering 
processes which lead the isolated Hubbard ladder toward the spin gap metal phase. Consequently 
when t± sets in, there exists, for any finite intraladder Hubbard repulsion, the region where the 
two-particle crossover dominates the one-particle crossover. 
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Appendix A: Derivation of the Scaling Equations 

The PRG program takes three steps: 1. coarse graining, 2. field rescaling and 3. renormalization. 
Coarse Graining 

As depicted in Fig. 18, first we divide the linearized intraladder bands with scaling-dependent 
bandwidth E(l) = E e~ l into the inner region, C um< (slow modes), and the outer shell, dC vm> 
(fast modes), where 

[k\\ || e vm {k\\) |< E(l + <H)/2}, 
dC vm> = {k {l | E(l + dl)/2 <\ e vm {k\\) |< £(/)/2} . 
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Coarse graining means integration of the fermion degrees of freedom over dC vm> with thickness 

1 



dE{l) |= -E(l)dl. 



(A-l) 



We leave e n and k± free. This procedure is performed for each branch (ym = LB, LA, RA, RB) 
and spin. 

The action (3.1) is then decomposed into the slow and the fast counterparts as 



S 



s < + s > 

(1) , c(2) 

< 



S 



S}'J+S^ + S ±< 



(2) ■ sf > + siS 



, c( 2 ) 



where S < consists of only the modes in C um< , while S > includes the modes in dC um> . The partition 
function is then rewritten as 



-/ 



T> < exp 



/ 



V > exp 



where T> < and P> symbolize the measure over the modes in C um< and dC um> , respectively. Under 



the weak coupling condition, t±,g^ <C t,t', we expand exp + + Sj^> + S 
regard to S^S i Sjo. ; an d 5*j_ 2 > > an d integrate out the modes in c?C> . Then we obtain 



?(2) 



7 (1) 



7 (2) 



with 



Z = Jv< exp 



slS + + 5« + + 



i,m,n=0 



where 



Iran 



l\m\n\ 



7 (2) 



7 (2) 



(A-2) 



(((' ' '))) = / ^> e '^O ' ') an d the subscript 'c' represents the connected diagram. Now we obtain 
the modified action, S < , with the modified intraladder propagators and interladder one- and two- 
particle hopping amplitudes written as 

£i(*n) = [i + *n(0<fl + o(di 2 )]g^(K), (a-3) 

t ± = [1 + 0(dZ 2 )]t ± , (A-4) 
gf = \1 + 4l(t)dl + 0(df)\gf , (A-5) 

^ M = + [*5,(0<« + ^ 2 )]^f> (A-6) 



where / symbolizes the generator of which is diagramatically shown in Fig. 7. In the low 
order arguments, z^ has [V^] n contributions with n = and n = 1. 

We obtain the 2^(0 > fff^P) an d -2^(0' by picking up and evaluating the Feynmann 

diagrams which give contributions in proportion to dl. The diagrams in Figs. 4(a) and 4(b) give 



(2) 

(/) and z^{l), respectively, which come from T2oo + ^00- The diagrams in Fig.4(c) give Z\\(l) 
which come from r 20 o- The diagrams in Fig.7 give f^f which come from r 2 2o- The second and 
third diagrams on the r.h.s of Fig.8 give Zjf„ which come from r 10 i + r 02- 
Field Reseating 

The right-moving sector of the coarse-grained intraladder kinetic action is 



EE 



E(l+dl)/2v F dk\\ 

-E(i+di)/2v F 2n 



1 + z^dl^XK^Rl^RUK). 



(A-7) 
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To restore the original cutoff and leave the intraladder kinetic action, Su, invariant, we tune the 
momenta, frequencies and field variables as 

K\\ = e dl K\\ = (1 + dl) K h 
R ma {K) = [1 + \{z\\{l) - 3}dl]R mcT {K), 

where 

K\\ = (e dl k h ie dl e n ), 
K = (e dl k b k ± ,ie dl e n ). 

We must leave the transverse momentum, k±, unchanged since the scaling law is well defined only 
along the longitudinal direction. Then we see that (A- 7) becomes 



r E{l)/2v F rfU,, 

EE / w ^^RumuK), 



-E(l)/2v F 2lT 

which is just the same form as the initial one, namely, scale- invariant. 
Renormalization 

Finally we renormalize the physical quantities in the action and set up the differential equations 
for them with keeping S\\ scale- invariant. Discarding all the 0(dl 2 ) terms in (A-3)-(A-6) and taking 
account of (A-6) and (A-8), we obtain the scaling equations 

^f=|(l)-2,a (A-8) 
dlnt ± (l) 



dl 



1-zM), (A-9) 



dV M (J) 

-^ = /^(0 + <(0^(0. ( A "!0) 

which give (3.9), (3.25) and Q3.29Q , respectively. Note that zf has [V^] n contributions with n = 
and n = 1. 

Appendix B: Evaluation of Self-Energy Diagrams 

We here give in detail the derivation of the factor, zn , which governs the one-dimensional crossover. 
In Ref.[13], the same results are given. The factor, Zn, comes from the renormalization of the 
intraladder propagotor, Q um (K\\). [y = R, L : m = A,B). Since Zn is independent of v and m, 
we consider only the case of Qrb\K\\). 

The 2nd order renomalization of Q um {K\\) originates from r 2 oo = \{{ S'm> ))• in Fig. 19, we 
give diagrammatic representations for I^oo which renoremalize Grb(K\\)- Although there are a lot 

(2) 

of ways to assign the outer-shell mode in On , there are only two possible ways which renormalize 
Qum{K\\). We specify these two cases by arrows in Fig. 19: 

case — I : kn G dCL mi> , k» + q\\ G Ci m2 <, k\\ — q\\ G Cij m3 <, 
case-E : k\ G C Rmi< , fcy + q\\ G dC Lm2> , k\\ - q\\ G C Rm3< . 

Then, among a lot of terms included in r 2 ocb the contribution to the renormalization of Qrb^v) 
is obtained as 



= 2 (2itv f ) 



^ 1)2 + 9 f-^9^)I BBB (dl) 
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+gf )2 X BAA {dl) + gfl AAB (dl) - gfgfl AAB {dl) 
+ gl 1)2 T AAB (dl) + gf )2 Z BA \dl) - g?g?X BAA {dl) 



where 



mim2m3 



(dl)=ir im2m3 (dl)+TZ lim2m3 (dl) 



with J™™ 3 (d/) and2^ im2m3 (rf/) corresponding to the cases I and E, respectively. Z™™ 3 ^/) 
is evaluated as 



dk[ f dq 

tdC Lrni> „ ^ 

de' Lmi [ dq T(e' Lmi ,k,q;T) 



=-[ 



'dc Lmi> 2irv F J 2-k D + 2v F (q - Ak F ) ' 
where Ak F = k Fmi — k Fni2 and D = ie n — v F (k — k F s)- The thermal factor is given by 
T(e Lmi ,K,q;T) = 



tanh £Lmi 



tanh 



£imi + v F (Ak F - q) 



2T 2T 
In the low-temperature limit, 

T(e Lmi ,k,q;T=0) = 



co th Vp ( q - A ^ - tanh Vp ^ k ~ ^ ~ 9 + 



2T 



2T 



4 for e Lmi [£ imi + v F (Ak F - q)\ < 0, 



for e Lmi [e Lmi + v F (Ak F - q)] > 0. 
Taking account of the restriction on available q\\ , coming from the case-I condition, we obtain 

1 



lp m2m:i (dl) = 



4tc 2 v f 



de' Tm . 17 dq+ I p de' Lrni 



E -\dE Q (l)\ Lrm y_B IL+Afc 



dq 



vp 



1 



D + 2v F (q - Ak F ) 

\dE (l) \ 1 JE -D][2E + D] 



16n 2 v 2 F [E + D] [2E - D] 
By taking the limit D / Eq <C 1 , we obtain 



X^ im2m3 (dl) 



D 



lQn 2 t 



-dl 



lQn 2 v 2 F 



<3RB{K)dl. 



Evaluation of X^ lim2m ' i (dl) gives the same result. Thus we obtain 
and the modified intraladder one-particle action 

K m=A,B o- 

where z\\ is given by (3.10). This result leads to (A-3). 
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Appendix C: Evaluation of Generators for Inter ladder Two-Particle Hopping 
Processes 

We here give a full account of evaluation of the generators for the two-particle hopping processes, 
f^f . It is convenient to rewrite the intraladder interaction part in terms of the composite particle 
variables as 



_ 71 U F v-^ LCDW /nBB* /nBB , CDW /nBA* /nAB , CDW /nBA* /n 
«Jint — 2-^i [9o ^CDW^CDW 9t ^CDW^CDW 9f W CDW W ' 

Q\\ 



i r ,SBW /nBB* /nBB _| ^SDWz-o-BA* /nAB , SDW/-0BA* ^BA (T 1 1 ^ 

"+" yo ^SDW ' ^SDW "+" 9t ^SDW ' ^SDW "+" 9f U SBW ' ^SDW \^' L ) 

i „SS/nBB*/nBB _i_ „SS ^>4 A _i_ ^SS/-/-)- 4 - 8 *^' 8 

~r yo ^ss ^ss + 9t ^ss ^ss + 9f u ss u ss 

, TS/^BB* AyBB _L „TS/*BB* (AAA , TS/nAB* /nAB} , / /i m 

where (9^ im2 denotes C^ im2 (Q||). 

As an illustration, we derive the generator for the CDW channel, /^ DW . In the action, (C-2), 
the CDW channel consists of 

CCDW _ 1 CDW /nBB* (n \/nBB (n \ 

^int - 2 Vpg 1^ U CDWW\\) LJ CDWW\\) 

Q\\ 

+ J"fE [^ DW ^cdwW||)0cdwW||) +^ CDW 0cdw(Q||)0cdwW!|)] + - B). 
Z Q\\ 

The generators depicted in Fig. 7(a) come from T 22- Among a lot of terms included in r 2 2o> the 
contribution to /^ DW , T^^, is obtained by assigning composite field variables to the outer shell 
mode as 

oCDW _ 1_„, „CDW r (/nBB* /nBB , /nBB* /nBB < /nAA* /nAA , /nAA* /nAA \ 
D int> — 2 Vp 90 Z"Q\\ l L/ CDW> L/ CDW< "+" L/ CDW< L/ CDW> "+" L/ CDW> Ly CDW< "+" L/ CDW< L/ CDW> J 

, ^ „, r„CDW//'/-) B - 4 * ^O- 4 ^ _l (nBA* /nAB , /nAB* /nBA , /nAB* /nBA \ 

-r 2 F ^Q\\ \-9t l^CDW>^CDW< ^CDW<^CDW> ^CDW>^CDW< ^CDW<^CDW>J 

, CDW / /nBA* /nBA , /nBA* /nBA , /nAB* /nAB , /-^AB* M 
"T 5/ l Ly CDW> Ly CDW< "T Ly CDW< Ly CDW> t l/ CDW> l/ CDW< + Ly CDW< Ly CDW>JJ 

+ (A<->fl), 

where 

0cDw< = r 1/2 E E EE^(*ii + Qii)W*ii), 

^cDw> = r 1/2 E E EEK^ii + Qii)W^ii)- 

Then we obtain 

pCDW _ 1//C-2 jqCDW\2\\ 
■■-220 — ^\W±>l D int> J // 



+ ^E [te CDW ) 2 ((^>^w>^w>)) + (^ CDW ) 2 ((5i>^w>^w>)) 

2 2 

CDW CDW V 1 / / q2 /nAB* (nAB \ \ /nBA* /nAB , / a m 
^ "T - y* y/ W±> Ly CDW> Ly CDW>// L/ CDW< Ly CDW< + l^ 1 ^ -°J> 

4 Q 



/„CDWn2 V^//c2 /nBB* r\BB \\/nBB* /nBB 
\9o ) W±> Ly CDW> ty CDW>// L/ CDW< Ly CDW< 



/nBA* /nBA 
L/ CDW< L/ CDW< 
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where the first, second and third lines on the r.h.s of the second equation correspond to the first, 
second and third lines of Fig. 7(a). We obtain, in the low temperature limit, 

(<OcDw>Si(9cDw>>> = -tlcos(q ± )T £ Y,G 2 Lm (h - 2k Fm ,ie n )G 2 Rm (k h ie n ) 

fc||Gc(CiJ m > £n 



2ttvf 
t 2 , dl 



t 2 ^ cos(q±d±) 



/ go(0 



E Q (l)-\dE Q (l)\ 
2 



+ 



go(0 

2 



E Q (l)-\dE (l)\ 
2 



de 



Rm 



fcoth" 2 ^ 2tanh%^ 

1 1 _|_ A 



'-2e RB y 



E (l)cosq±. 



TTVf 

Similar manipulation gives 



\\ Ly DW°± Ly DW // 



t\dl 

71 Vp 



E 2 (I) cos q±, 



where m and m denote different bands. Thus we obtain 

F CDW _ tjTVVp . CDW\2 \- sr\BB* n BB 

1 220 — ~ 8E 2 (l) COS y- LL/ CDW< Ly CDW< 



+ 8^(Z) 



(ft CDW ) 2 + (% CDW ) 2 ]Ecosgi«<^ 



CDW< 



^±/_£_ jj CDW CDW /r\BA* n AB 

4£/ ol<J Q,, 



After the field rescaling procedure, the first, second and third lines of the above equation give 
the first terms in (3.30), (3.32) and (3.33), respectively. It is here useful to present the following 
results. 



/ / f~) mm Q2 f)inm* 



}} = - t £E^ 2 (l)cosq ± , 
((OfiftSlOSF)) = +^£ - 2 (Z)cosg ± , 
((O^S 2 ± 0^*)) = +^£- 2 (/)cosg ± , 
((O^SiOr*)) = -^E 2 (l)cosq ± , 

where DW=CDW/SDW and S=SS/TS. When we evaluate the second and third diagrams of 
Figs. 8(a) and 8(b), it is useful to note the outer shell integration, which appears in the evalua- 
tion of Ton and roo2> is given by 



((o M 

for M = CDW, SDW, SS or TS. 



mm/r\mm*\\ / / //ymm //ymm* 

U M II — \\ U M u 



'M 



dl 
TTVf 
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Fig. 1: Array of Hubbard ladders studied here. 

Fig. 2: Schematic illustrations of the one-particle process and the two-particle process (in the d- 
wave superconductivity channel). In the one-particle process, a particle hops from one ladder to 
a neighboring one, while in the two-particle process, a pair of particles (bipolaron) hops from one 
ladder to a neighboring one. 

Fig. 3: (a) Four branches (Li? , LA, RA, RB) of linearized bands with the bandwidth cutoff E, 
and (b) intraladder two-particle scattering vertices The solid and broken lines represent the 
propagators for the right-moving and left-moving electrons, respectively, m and fh denote different 
bands. 

Fig. 4: Diagrams which contribute to the 3rd order vertex corrections, (a) to the intraladder 
backward scattering processes g^}\ (b) to the intraladder forward scattering processes gj?\ and 
(c) to the intraladder self-energy processes. Wavy lines represent one of the intraladder scattering 
vertices, gf (i = 1, 2; fi = 0, /, t). 

Fig. 5: Diagrammatic representation of the scaling equation for the inter ladder one-particle hopping 
amplitude. A zigzag line represents the interladder one-particle hopping, t±. 

Fig. 6: Scaling flows of the interladder one-particle hopping amplitude, t±, (a) for U = 0.3, and 
(b) for t_i_o = 0.01. For smaller U and larger t±o, t±(l) exceeds unity in the course of the scaling. 

Fig. 7: Diagrammatic representations of the generators of the interladder two-particle hopping 
processes, fjjf, (a) for the density wave channels, and (b) for the superconducting channels. The 
first, second and third lines correspond to the flavor indices ji — 0, / and t, respectively. Hatched 
circles represent the coupling strengths, g^f . 

Fig. 8: Diagrammatic representations of the leading order scaling equations for the two-particle 
hopping amplitudes, V^ 1 , (a) for the density wave channels (M =CDW or SDW), and (b) for the 
superconducting channels (M =SS or TS). Hatched squares represent V^f . 

Fig. 9: (a), (b), (c) Scaling flows of the intraladder coupling strengths for the composite fields, 
g^f , and (d), (e), (f) scaling flows of the interladder two-particle hopping amplitudes, Vj^ , for 

U = 0.2, 0.3, 0.4, respectively, and t±o = 0.01. The vertical broken line corresponds to the 
scaling parameter, l c , at which V SCd diverges. 

Fig. 10: Scaling flows of the stiffness of the total-spin mode, K a+ , the interladder one-particle 
hopping amplitude, t±, and the interladder two-particle hopping amplitude in the SCd channel, 
V SCd , for various initial conditions (U,t±o). For the region, I > Min(Z c , / cross ), the scaling flows, 
drawn by broken curves, have no physical meaning, since the weak coupling picture breaks down 
in the region. 

Fig. 11: Phase diagram of the weakly coupled Hubbard ladder system spanned by t±o and the 
reduced temperature T = T/E for U = 0.3. SGM, SCd and 2D denote the sping gap metal 
phase, the d-w&ve superconducting phase and the two-dimensional phase, respectively. Gradual 
chanage of darkness in the SGM phase schematically depicts the gradual approach of the isolated 
systems to their low-energy asymptotics. Thick broken lines denote the crossover boundaries. 
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Fig. 12: Phase diagram of the weakly coupled Hubbard ladder system spanned by the intraladder 
Hubbard repulsion U and the reduced temperature T = T/Eq for tj_o = 0.01. Gradual chanage 
of darkness in the SGM phase schematically depicts the gradual approach of the isolated systems 
to their low-energy asymptotics. Thick broken lines denote the crossover boundaries. 

Fig. 13: Dependence of the crossover value of the interladder one-particle hopping amplitude, t± c , 
on the intraladder Hubbard repulsion, U. For tj_ > t± c , the system undergoes the crossover to 
the 2D phase, while for tj_o < t± c , the SGM phase transits to the SCd phase. 

Fig. 14: (a) Coupled ladder system with the misfit of the rungs in the neighboring ladders 
which actually exists in real Sri4_. E Ca :E Cu2404i compounds and (b) static interband polarization, 
Xab(q\\,Q±), on the 1st Brillouin zone along the line T(0,0) — > X(7r, 0) — > M(7T,7r) — > T for 
t"/t = 0, 0.2, 0.3 and temperature, T = 10~ 5 t. 

Fig. 15: (a) Two branches(L, R) of linearized bands with the bandwidth cutoff E, and (b) intrachain 
two-particle scattering vertices The solid and broken lines represent the propagators for the 
right-moving and left-moving electrons, respectively. 

Fig. 16: Scaling flows of the interchain one-particle hopping amplitude, £ c hain± (the upper half 
plane ), and the interchain two-particle hopping amplitude in the SDW channel, \/ SDW (the lower 
half plane ), for U = 0.1, 0.2, 0.3, and t ±0 = 0.01. 

Fig. 17: Phase diagrams of the weakly coupled chains, (a) for U = 0.3, and (b) for t c hain±o = 
0.01. TL and 2D denote the Tomonaga-Luttinger liquid phase and the two-dimensional phase, 
respectively. Thick broken lines denote the crossover boundaries. 

Fig. 18: Division of the four branches of the linearized bands into the inner region, C vm< , and the 
outer shell, dC vm> . 

Fig. 19: All labeled diagrams which renormalize Grb(K\\)- Possible combinations of the band 
indices are m 1 m2ro3 = BBB, ABA, BBA. Inserted arrows with the labels I and H denote the 
outer-shell modess in the cases I and E, respectively. 



25 



t 









t' 



















Fig. 1 



One-Particle Process Two-Particle Process 




Fig. 2 




Fig. 4 




2 4 6 8 



Fig. 6 



(a) CDW or SDW 



O-process 




m/O v ^m/0\ 



m 







<^ m \ y 






f-process 








/ s 



m /f\m/f V m m^^VH^^ 



m 



t-process 




m 



f ^ynyt\ 



m m\ 



m m 




t\m/f \- 



m 



(b) SS or TS 



O-process 








m \, 


f-process 







t-process 




m m\ 



m 



m 




SDWf 
SDWt 



/V/ 

T 

.02 



01 - 



T * 

cross 



SGM 




T 



c 



SCd 



2D 



0.01 0.02 f 0.03 

txc 

Fig. 11 



0.04 

t-LO 



-2 



xiO 
1 

. 8 
. 6 
0.4 
0.2 



r t 



cross 



2D 



SGM 




SCd 



T 



C 




0.1 0.^ 0.3 0.4 

Uc 

Fig. 12 



. 5 



(a) (b) 

L R R 




R L R 

Fig. 15 



0=0 .1,0.2, 0.3 




U=0 .3,0.2,0.1 



Fig. 16 




0.01 0.02 0.03 0.04 



t chain _L0 




T T 



LB LA 



RB 

ZpCy m> 




Cym< 



Vm> 



Fig. 18 



Lmi,K f 



i 



JL 



Lm2,K f +Q, 



RB, K 



2 



RB, K 



Rm3,K-Q 



Lmi,K f _ 




Rm3, K-Q 

RB,Kl..^i <?RB, K 

Lm2,K f +Q 




Lm2,K'+Q Lmi, K ' Rm3, K-Q Rm3, K-Q Lmi, K ' Lm2,K'+Q 



Fig. 19 



